Density of states of the interacting two-dimensional electron gas 
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We study the influence of electron-electron interactions on the density of states (DOS) of clean 
2D electron gas. We confirm the linear cusp in the DOS around the Fermi level, which was obtained 
previously. The cusp crosses over to a pure logarithmic dependence further away from the Fermi 
surface. 

PACS numbers: 



It was established more than 20 years ago by Altshuler 
and Aronov 1] and Altshuler, Aronov, and Lee that 
in low-dimensional diffusive systems the electron-electron 
interaction leads to suppression of the single-particle den- 
sity of states (DOS) at the Fermi level. In two dimen- 
sions their theory predicted logarithmic cusp at the Fermi 
level due to the diffuson pole divergences in the vertex 
renormalizaion. The theory was later extended by Rudin, 
Aleiner and Glazman |3( to larger values of the energy 
measured from the Fermi surface. 

Recently it was shown by Khveshchenko and Reizer |4| 
and by Mishchenko and Andreev [5( using the diagram- 
matic approach with the RPA dynamical susceptibility, 
that large electron-electron interaction induced correc- 
tion to the DOS exists even in the absence of disorder. 
Both groups obtained a linear cusp at the Fermi level 
8v(e)/vn ~ |e| / Ep (with the slope differing by a factor of 
2), independent of the strength of the electron - electron 
interaction r s . More recently Rollbuhler and Grabert 
using a path integral technique (applied for relatively 
small couplings ?' s < 1) found numerically that the slope 
does depend on the coupling and flattens away from the 
Fermi surface. 

Due to apparent discrepancy between the two ap- 
proaches we reconsider the problem in the framework of 
the diagrammatic method. Our aim is to spell explicitly 
and to clarify all the approximations which are made to 
obtain the transparent analytical formulas for the DOS 
correction. We also obtain the results which are partially 
different from those of Refs. 0, 0. 

The Hamiltonian of the 2D electron gas is 



H = Y1 £ P a U a P° + 9 a l+q<7 a i>- q <j> v {<l) a p'<y' a P< T i 
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The Green's function G(p, E) is given by the equation 
G- 1 (p,E) = Go\p,E)-X(p,E), (3) 



where Go(p,E) is Green's function in the absence of 
electron-electron interaction 



G^E) = E-L-; 



(4) 



the self-energy S(p, E) in this paper will be calculated in 
the Random Phase approximation (RPA) at T = 0: 

^> = /S/£c„(p + ,£ + .)ia (5) 

The DOS is 

v{E) = -~lmJ^G(p,E + iQ). (6) 

Expanding Green's function with respect to the self en- 
ergy and reversing the order of integrations, one obtains 
the electron-electron interaction induced correction to 
the DOS Q 

8v{E) = ~Im J ^2-SG(p, E + i0) = ImX(E + iO). 

(7) 

We will work in Matsubara formalism, calculate X for 
imaginary frequency 

d 2 q f°° dej V(q) 
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and at the end will make an analytical continuation. 

Using the fact that Gg(p, ifl) = idGo(p, ifl)/dtt, we 
can write the integral over p of the three Green's func- 
tions as 



d2 P r-,-2 



(2tt)< 



G' (p,in)G Q (p + q,i(n + Lu)) 
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where the polarization function 11(0, q, uj) is given by the 
equation 



II(g, iw, ifl) 
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The last equation in polar coordinates takes the form 
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Since the main contributions come from fermionic mo- 
menta close to Fermi momentum, one can make the first 
approximation by replacing p in the term pq cos((9) / 'm by 
\/2m/i = mvp and ignoring the term q 2 /2m. After that 
the integration over p can be easily performed: 



H(q, iw, ifl) 



dO 



(2tt) 2 Jq iw — VFqcos(9) 
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where 



L(w, Q., q) = log 



(fi + w) — iji + ivFqcos{9) 



Q — ifj, 

Explicitly, presenting the logarithm as 

1 (n + w) 2 + ^-v F qcos(9)) 2 
L(w, 0, «) = - log ^—^ 



(13) 
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and taking into account that f2, w, vpq <C /i, we make the 
second approximation: 

L = ni [0(-fi)8(fi + w) - 8(0)8(-fi - w)} . (15) 

After that, the integral over 9 in Eq. (|12fl can be easily 
calculated, and the polarization function takes a form 



H(q, iw, ifl) 



m sign(w) 



2 y/w 2 + v^q 2 
[9(-n)6(fi + w) - 8(fi)8(-Q - w)] . 

Taking appropriate derivatives one obtains 



(16) 
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exact Eq. I|10(l and (as it is traditionally done) integrate 
over f2 first, to obtain 



P(q,iu;) = 2j2- 
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where n p is the Fermi distribution function. Obvious 



algebra gives 
P(q,iw) -- 
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Integrating over p, 
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and expanding the radical in small g, one obtains 
2e 2 m 
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y/w 2 + v' F q 2 
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Now we are ready to return to the calculation of the 
DOS. Substituting Eqs.|T7|) and Q into Eq.©, consid- 
ering for defmiteness > 0, and subtracting an "inessen- 
tial" constant, one has 



X(iQ) 



2tt 2 



qdq 
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(24) 



W ,-uvr,-/-"vr ■ *>-v- ■ -„ 2 [^ 2 + i; 2 ,g 2 ]3/ 2 

[e(-n)e(n + w)-e(n)e(-n-«)]. (17) 

The dielectric constant e(q,w) is related to the polar- 
ization operator 



Apart from the fact that we are considering imaginary 
energy, Eq. i|24|) coincides with those obtained in Refs. 
0, 0. (Comparing with the Eq. (2) from the latter ref- 
erence, we have an additional factor 1/2.) We notice, 
however, that the double integral can be calculated in a 
more rigorous way than it was done there. After we intro- 
duce the dimensionless variable q = vpqjw, and change 
the order of integrations, the integral takes a form 



P(q, iw) = -2 J |-n(g, iw.ifl) 



by the equation 

e(q,iw) = l-V(q)P(q,iw). 



(18) 
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It is easier, however, not to use the approximate Eq. illtj|) 
for the polarization function, but to insert in Eq. 1181) 



The integral over w can be easily calculated and we obtain 

xm e2 ™ r d<} 



2ir 2 v F J Q [1 + q 2 f/ 2 
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To obtain the density of states we should substitute — > 
— ie, where e = E — \i is the energy measured from the 
Fermi surface, and take the imaginary part. Thus we 
obtain 



dq 



Atzvf Jo [i + <f] 3/2 



log f(e,q), (27) 



where 



f(e,q) = 1 



2 -9 

x q 



(28) 



and 



2me 2 VF ' 



(29) 



and — m/ir is the DOS of the noninteracting 2DEG. 
For i> 1 we may put in Eq. (|27|) 



f(x) = 21ogx. 



Hence for |e| > (e 2 /v F )E F 
6v(e) e 2 



log 



2ttvf 2me vf 



(30) 



(31) 



For i«l the main contribution to integral in Eq. H27|) 
comes from the vicinity of the lower limit, hence we can 
expand all the functions with respect to q, thus obtaining 
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The integral can be easily calculated and we obtain for 
|e| « (e 2 /v F )E F 



Me) 



SEp 



(33) 



The equations (|31|) and (|33|l are our main results. The 
correction is smaller than v§ consistently with the per- 
turbative assumption. At small |e| our result coincide 
with that of Ref. up to a factor 1/2 (see the state- 
ment immediately after Eq. (|24(l 'l. The linear segment 
crosses over to the logarithmic one at the energy scale 
|e| = (e 2 /vf)Ef- This last statement is in good agree- 
ment with numerical results of Ref. 0] . 

We are grateful to E. G. Mishchenko and M. Reizer for 
their criticism which helped us to find a mistake in our 
previous version. 
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